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Charge and spin transport in a junction involving two triplet superconductors and a ferromagnetic 
barrier are studied. We use Bogoliubov-de Gennes wavefunctions to construct the Green's function, 
from which we obtain the Josephson currents in terms of the Andreev reflection coefficients. We 
focus on the consequences of a finite barrier width for the occurrence of O-tt transitions and for the 
spin currents, and examine the appropriateness of the common ^-function approximation for the 
tunneling region. 
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I. INTRODUCTION 

Josephson junctions between spin-singlet superconduc- 
tors with ferromagnetic tunneling barriers have recently 
attracted much interest. Since spin-singlet supercon- 
ductivity is strongly suppressed by the exchange split- 
ting of the Fermi surface in the ferromagnet, such junc- 
tions represent a unique way of probing the interplay of 
two competing states. A startling consequence is the 
so-called O-tt transition: By varying the width of the fer- 
romagnetic tunneling barrier or changing the magnitude 
of the exchange splitting, the current vs. phase relation- 
ship reverses sign relative to the non-magnetic case.^^^ 
Roughly speaking, the pair-breaking effect of the ex- 
change fields causes a spatial oscillation of the super- 
conducting pairing correlations within the ferromagnetic 
tunneling barrier, which can therefore connect the two su- 
perconductors with either a or tt phase shift. This phe- 
nomenon has been the subject of intense theoretical^^^ 
and experimental^^^^ interest and is by now very well 
understood. 

At least from a theoretical point of view, it is there- 
fore surprising that the physics of triplet superconductor 
(TSC)-ferromagnet (FM) bilayers^^^^ and junctions has 
received relatively little attention. ^^'^^'^^^^ Due to the 
intimate relationship between triplet pairing and ferro- 
magnetism, one would anticipate very different behavior 
compared to the singlet case. Recently, several papers 
have pointed out that the orientation of the magnetiza- 
tion relative to the TSC vector order parameter is a key 
parameter of such systems, leading to a variety of exotic 
effects: A universal O-tt transition dependent upon the 
orientation of the magnetic moment, an intimate connec- 
tion between spin and charge Josephson effects, and also 
phase- independent Josephson spin currents. ^^'^^'^^^^ Al- 
though none of these predictions have yet been verified 
experimentally, the recent creation of superconducting 
thin films of Sr2Ru04 is a significant step towards the 
fabrication of TSC-FM heterostructures.^^ A deeper un- 
derstanding of the physics of these systems is therefore 
desirable. 

Most investigations of TSC junctions were performed 
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FIG. 1: Schematic diagram of the junction. 



assuming an atomically thin barrier modeled by a ^- 
function separating the superconducting regions. The 
effect of a finite barrier width have only been studied by 
Rahnavard et al.^'^ for the case of a {pz + ipy)-wdive gap 
where the authors obtained results that are only partially 
consistent with the previous perturbative calculations^^ 
for a ^-function barrier. This raises the interesting possi- 
bility that additional mechanisms for the Josephson cur- 
rents appear due to the finite barrier width. It is there- 
fore the aim of this work to systematically investigate the 
relationship between the results in the (5-function limit 
and for a more realistic barrier model. In our calculations 
we utilize a non-perturbative method to determine the 
Josephson current by solving the Bogoliubov-de Gennes 
(BdG) equation. We compare the results with the pre- 
dictions of perturbation theory.^^ 



II. THEORY 



We consider a Josephson junction consisting of an FM 
layer of width L sandwiched between two TSCs, shown 
schematically in Fig. 1. The BdG equation for quasipar- 
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tide states with energy E reads 
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The noninteracting Hamiltonian is 
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where for simphcity we assume two-dimensional circular 
Fermi surfaces lying in the y-z plane, as well as equal 
effective masses and chemical potentials in all three re- 
gions of the junction. The magnetic moment M of the 
FM lies in the x-y plane and forms an angle a with the 
X axis so that 



M = M (cos a, sin a, 0) . 



(3) 



This formulation is completely general since the TSCs 
are invariant under spin rotations around the x axis. 
The gap matrix in equation (1) is 



A{r) = ALe{-z) + AnOiz - L), 



(4) 



where {v = L, R) is related to the d vector by A^ = 
i&y (<T • d^). We assume the TSCs to be in an equal-spin- 
pairing unitary state. In particular, we take 
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A(k)e, 
A(k)e^' 



(5) 
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where A(k) denotes the orbital pairing state, so that the 
two superconductors differ only in the overall phase (j) of 
the order parameter. In this work we will be concerned 
with the three orbital pairing states 



A(k) 



A{T)ky/kF, 

A{T)k,/kF, 

A{T){k,^iky)/kF, 



p^-wave, 

p^-wave, (7) 



where A(T) ^ Aq^{Tc -T)/Tc is the weak-coupling 
temperature-dependent gap magnitude and kp is the 
Fermi momentum of the normal state. 

The BdG equation is solved separately in the differ- 
ent regions under the assumption of low-energy excita- 
tions, i.e., /i |A(k)|, so that the electron and hole 
wavevectors can be regarded as approximately equal in 
magnitude, ke ^ kh ^ kp-, the so-called Andreev approx- 
imation. The dephasing of the electron and hole wave- 
functions in the finite FM layer constrains the validity 
of this approximation to {ke — kh)L ^ 2EL/hvF <C 1, 
where vf ^ lO^ms"^ is the Fermi velocity and E < 
max{|Ak|} ^ 0.1 meV (we assume Tc ^ IK). This 
implies that we restrict ourselves to thin layers with 
Lkp ^ 100. In the TSC we obtain plane-wave solu- 
tions '^k,e(/i),cr = <f>\i,eih),a^^^'^ electrou-like (hole-like) 



quasiparticles with wavevector k and spin cr. The spinors 
are 



0k,e,t = (5tk^k,0, Vk,0)^, 
0k,/i,t = (5tk^k,0,2ik,0)^, 
0k,e,i = (0,5;k^k,0, Vk)^, 
^k,/i,i = (0,S;k^k,0,iik)^, 



(8) 
(9) 
(10) 
(11) 



where ^^k = y^{ETn^/2E, ^k = ^/{E - ni^)/2E, 
^^ k = V^ ^-|A(k)P, 8,k = -aA(k)/|A(k)|, and |k| = 
^J2m|l/h? = kp. In the FM region we obtain plane- wave 



= ^f^^ s^^^^ ^ electron-like (hole- 



solutions ^Is^^^y,^^, 

like) quasiparticles with wavevector k^. The spinors of 
spin 5 = ± electrons and holes are 



<f = («e-'7V2,l/V2,0,0f, 
0f,f = (0,0,se-/\/2,l/y2f. 



(12) 
(13) 

Due to the exchange splitting in the FM, we have 
majority-spin and minority-spin Fermi surfaces with radii 
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(/i + sqiibM) = krVlTsX, (14) 



where A denotes the ratio of the exchange splitting to the 
chemical potential. 

We use these piecewise solutions to construct the wave- 
function ansatz for a spin-a electron-like quasiparticle in- 
cident from the left TSC with wavevector k = {ky^ kz) = 
/cF(sin6>, cos6>), where 6 is the injection angle relative to 
the z axis: 

*e..(r) = 

cr'=n 

s=+- 

+ e{z - L){ J2 [daa'^t,^' + ^-'V^-k,.,.']}- (15) 

A schematic representation of the wavevectors and trans- 
mission and reflection amplitudes appearing in the wave- 
function ^e,cr(r) is shown in Fig. 2. For z < this ansatz 
describes Andreev reflection of hole-like quasiparticles 
with wavevector k and specular reflection of electron- 
like quasiparticles with wavevector — k = (ky^—kz). 
These processes are weighted by their probability am- 
plitudes a^^/ and b^^,^ respectively. In the FM (i.e., 
for < z < L), the wavevectors of the transmitted 
and reflected quasiparticles are spin-dependent. For 
wavevectors of spin-s right-moving (left-moving) elec- 
trons (holes) we write k^ = (/cy,/c|), while spin-s left- 
moving (right-moving) electrons (holes) have wavevec- 
tors — k^ = {ky^—kl). The associated probability am- 
plitudes for electrons (holes) are c^^ and e^^ (d^^ and 
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FIG. 2: (Color online) Schematic representation of the various wave vectors, amplitudes, and coefficients for a spin-cr electron-like 
quasiparticle incident from the left TSC, with wavefunction ^e,cr(r). The arrows denote the directions of travel of electron- like 
quasiparticles (solid blue) and hole- like quasiparticles (dashed red). 



f^s)- -^o^ z > L we have transmitted electron-like and 
hole-like quasiparticles with wavevectors k and — k, and 
probability amplitudes g^^, and h^^^^ respectively. The 
ansatz for an incident hole-like quasiparticle ^/^^cr(r) is 
analogous. The probability amplitudes for this case are 
distinguished by the superscript e.g., a^^,, b^^, etc. 
Due to the translational invar iance along the y axis, the 
wavevector component parallel to the interface is con- 
served during scattering, which gives 



sin^ = kp sin^^. 



(16) 



where 0^ is the transmission angle for quasiparticles from 
the spin-5 Fermi surfaces. Making use of this relation, k^ 
can be expressed as 



K = cos6>^ = /cfa/cos^ + sX. 



For injection angles 



> c^crit = arccos 



(17) 



(18) 



is imaginary for minority-spin quasiparticles, i.e., the 
wave decays exponentially into the FM region. 

The probability amplitudes for the various processes 
are calculated from the continuity of the wavefunction 
and its derivative at the interfaces. 



^e(/i),a(r)U=o- = ^e(/i),a(r)U=0+, 
^e(/i),a(r)U=L+0- = ^e(/i),a(r)U=L+0+, 



^^e(.),.(r) = ^^e(.),.(r) 



z=L-\-0- 



^^e(/^),a(r) 



z=0+ 



(19) 
(20) 

(21) 

(22) 



The results for a (5-barrier can be obtained by taking the 
limits L ^ 0, A ^ OO while keeping the product LA 
constant. 

The charge current Ic and spin current (/i = x^y^z) 
perpendicular to the junction interface are obtained using 
the Furusaki-Tsukada formulas. ^^^^ These are written in 
terms of the Andreev reflection coefficients a^^/ and a^^, , 
where the energy argument is analytically continued to 
Matsubara frequencies [E iuon = i (2n — l)7r//3]. 
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X < 7^ <a(k,CJn) 



n, — k 



at(-k,u;„)^ (23) 



X J lM„e^(k,^„) _ l^„^^(_k,^„) \ , (24) 



n, — k 



Is = - j '^^yHH- 



X JlM«e_^(k,c.„)-l^a^,(-k,^„)l, (25) 



where fln,h = \/^n + l^kp, = — cr and / dky denotes 
the integration over all momenta parallel to the inter- 
face.^^ The X component of the spin current vanishes 
since the Cooper-pair spin is always perpendicular to the 
d- vector. 
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FIG. 3: Phase diagram of the junction showing 0- and vr-states as a function of the angle a and barrier width L for the three 
considered gap symmetries. In all panels we set T = 0.4 Tc and A = 0.3. 
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FIG. 4: Phase diagram of the junction showing 0- and vr-states as a function of the angle a and exchange splitting A for the 
three considered gap symmetries. In all panels we set T = 0.4 Tc and L = 20k . 
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FIG. 5: Phase diagram of the junction showing 0- and vr-states as a function of the exchange splitting A and the barrier width 
L for the three considered gap symmetries. In all panels we set T = 0.4 Tc and a — 0. 



III. RESULTS 

We have numerically solved Eqs. (19)-(22) for the An- 
dreev reflection coefficients at T = 0.4 T^. At least 400 
Matsubara frequencies have been used in all calculations. 



A. Charge current 

We are mainly interested in the occurrence of O-tt tran- 
sitions, which involve a change of the phase difference 
of the global free-energy minimum between ^ = and 



(j) = TT. Points in parameter space can be identified as 
corresponding to a 0-state or a 7r-state by examining the 
current- vs. -phase relationship. Specifically, we have the 
thermodynamic relation^^ 



(26) 



SO that /c = implies an extremum of the free energy 
F. In our system, minima only occur at ^ = or tt. We 
therefore numerically calculate the derivative of IcW at 
these two points to find out which one is a minimum. If 
both are minima we integrate IcW to obtain the free- 
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FIG. 6: Critical charge current as a function of the barrier width L. In all panels we set T = 0.4 Tc and A = 0.3. 




FIG. 7: Critical charge current as a function of the exchange splitting A. In all panels we set T = 0.4 Tc and L = 20 /c^^. 




FIG. 8: Critical charge current (black solid lines, left axis) and integrated transmissivities Tint (red dashed lines, right axis) as 
functions of the barrier width L. In all panels we set a = 0, T = 0.4 Tc, and A = 0.3. 



energy difference AT between = and tt, 

AF=^£ dcpic{4>), (27) 

which allows us to identify the global free-energy mini- 
mum. Using this method, we construct phase diagrams 
of the junction for cuts through the (a, A, L) parameter 
space, which are presented in Figs. 3-5. 

We find that O-tt transitions are possible for all consid- 
ered gap symmetries by increasing the magnitude of the 
magnetization, rotating the magnetization within the x-y 
plane, or by increasing the barrier width L. For A < 1 in 
the (a. A) maps and in the (a, L) maps, we find multiple 
dome-shaped regions where the junction is in a 7r-state. 
These 7r-state domes are only found for a < 7r/4 and their 
appearance is shifted to lower values of a by increasing 



L or A. This effect is strongest for the p^-wave case, 
while the tt states are much more robust for Pz-wa^e and 
{Pz + ipy)-w8iYe pairing. In the (a. A) maps the 7r-states 
are strongly suppressed for A > 1. 

For a = 0, the distribution of 0- and 7r-states in Fig. 
5 shows a stripe-like pattern, which does not strongly 
depend on the pairing symmetry of the TSC but rather 
upon the characteristics of the barrier. In particular, the 
sharp bend of the stripes at A = 1 in all panels is a 
signature of the onset of half-metalicity since |k^| van- 
ishes at A = 1. Interestingly, these phase diagrams are 
very similar to those obtained for a junction involving 
singlet superconductors on both sides of the barrier. -'^^ 
The parallels between these different junctions originate 
from the universal pair-breaking effect of the FM-layer. 
Note that for the TSC case the FM layer is purely pair- 
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breaking only for a = 0, and pair-breaking is completely 
suppressed upon increasing a to a = 7r/2. In contrast, for 
a singlet-superconductor junction the FM layer is purely 
pair-breaking for all a. A O-tt transition caused by vary- 
ing the orientation of the magnetic moment is therefore 
an unambiguous experimental signature of triplet super- 
conductivity. 

The direct observation of 7r-j unctions is not trivial, as it 
requires the detection of half vortices at the barrier. ^^'^^ 
A more practical experimental signature is provided by 
the critical charge current 



of the gap. We therefore define the integrated transmis- 
sivity 



;{|/c(0)|}, ^G[0, 



(28) 



As the junction is tuned through a O-tt transition, the sign 
change of the dominant sin (j) term in the charge current 
/cl^) causes the appearance of sharp, cusp-like minima 
where /crit nearly vanishes.^ If /c(0) were exactly pro- 
portional to sin0, this sign change would exactly coin- 
cide both with /crit = and with a O-tt transition; the 
presence of weak higher harmonics in /c(0) complicates 
the analysis, but the minima in /crit nevertheless remain 
a good proxy for the O-tt transition. We plot /crit as a 
function of L and A for a = 0, 7r/4, and 7r/2 in Figs. 6 
and 7. The sharp minima in the critical current charac- 
teristic of the O-TT transition are only found at a = 0, and 
are in good agreement with the phase boundaries in Figs. 
3 and 4. Upon increasing a, the sharp minima disappear 
by a = 7r/4 and are replaced by broad minima in /crit, 
which coincide with every second maximum of the a = 
critical current. These minima do not correspond to sign 
changes of the sin0 term. At a = 7r/2, most structure in 
/crit has vanished. 

The charge current found by perturbation theory for a 
(5-barrier junction is,^^ 



/g,^'*(0) oc (Tsp - Tsf cos2a) sine 



(29) 



where Tgp and Tgf denote the positive tunneling matrix 
elements for spin-preserving tunneling and spin-fiip tun- 
neling, respectively. As one can see, a 7r-state can only 
be realized if Tgp < Tgf and a < 7r/4, with the latter 
condition clearly consistent with the phase diagrams and 
critical current plots. 

It is interesting to see if Eq. (29) is also valid in the 
finite-barrier case if we regard Tgp and Tgf as functions 
of L and A. The matrix elements Tgp and Tgf are ex- 
pected to be roughly proportional to the corresponding 
transmissivities.^^ Note that the transmissivities are in- 
dependent of a for both the (5-function and finite width 
barriers. From Eq. (29) we see that the spin-fiip trans- 
missivity T^f should be larger than the spin-preserving 
transmissivity T^p in a 7r-state, whereas the opposite 
should be true in a 0-state. We have obtained T^f and 
T^p from the scattering matrix of the FM barrier, which 
does not change below Tc^'^^ As the Josephson current 
is expected to be proportional to the square of the TSC 
gap magnitude, we introduce an angular weighting of the 
transmissivities to account for the direction dependence 



with 



A, = 




for p^-symmetry, 
for p^-symmetry, 
for pz + ipy-symmetTy . 



(30) 



(31) 



Spin-fiip tunneling is dominant for Tint > 0.5, whereas for 
Tint < 0-5 spin-preserving tunneling is most important. 
As seen in Fig. 8, the barrier widths for which Tint = 
0.5 are in very good agreement with the sharp minima 
of the critical current at a = and thus with the 0- 
TT transitions. This shows that Eq. (29) gives a good 
description of the relevant physics if one includes the L- 
and A-dependence of the transmissivities. 

In Fig. 8 we also see that the convergence of Tint 
towards 0.5 with increasing L is much faster for py- 
symmetry compared to the other cases. Due to the Py- 
wave form factor, shallow incident angles dominate the 
integral in Eq. (30), which implies averaging over quasi- 
particle trajectories corresponding to a large range of ef- 
fective barrier widths. Since the probability of spin-fiip 
transmission oscillates with the effective barrier width, 
increasing L rapidly leads to equal integrated probabili- 
ties of spin-fiip and spin-preserving transmission. In con- 
trast. Tint decreases more slowly for p^-symmetry since 
the weighting favors normal incidence, and thus the av- 
erage is over a smaller range of effective barrier widths. 
The angle-independent weighting for the {pz -\- ipy)-wdive 
case lies between these two extremes. From Eq. (29) we 
hence deduce that the tt state in the Py-waye junction is 
much less robust to increasing L than for the other cases, 
consistent with Fig. 3. Since increasing A decreases the 
length scale of the oscillation between spin-fiip and spin- 
preserving transmission, this argument also explains the 
strong suppression of the tt states in the Py-waye junction 
at moderate to strong polarizations in the (a. A) maps 
(Fig. 4). 

The loss of every second peak in the a = critical cur- 
rent upon increasing a to 7r/4 can be understood from 
Eq. (29) as resulting from a vanishing spin-fiip contribu- 
tion. As evidenced by the plots of Tint in Fig. (8) the 
relative strength of T^f and T^p oscillates with L and A. 
By removing the spin-fiip contribution at a = tt/A, the 
variation of the spin-preserving transmissivity on these 
quantities is evidenced in the critical current. Further 
support for this interpretation comes from the observa- 
tion that every second maximum in the a = critical cur- 
rent also corresponds to a maximum in Tint, i'C, strong 
spin-fiip tunneling and weak spin-preserving tunneling. 
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FIG. 9: Spin-current components (a) /|, which is due to spin 
filtering, and (b) J^, due to spin flip, as functions of the phase 
difference 0. In both panels we set L = 10 kp^, A = 0.3, and 
T 0.4 Tc. 



B. Spin current 

We find that the polarization of the spin current al- 
ways lies in the y-z plane. The y and z components are 
plotted in Fig. 9 for the p^-wave case; the results are 
qualitatively identical for the other orbital symmetries. 
We show below that these two polarizations are due to 
different mechanisms, as proposed by Brydon et al.^^ 

The y component of the spin current shows the a and (j) 
dependence expected for a spin- filtering mechanism. ^^'^^ 
It occurs when the magnetic moment of the barrier has a 
component perpendicular to both d- vectors of the TSCs. 
The spin-t and spin-^ states with respect to this direc- 
tion are not mixed by scattering at the barrier but do 
experience different barrier transmissivities, leading to 
spin filtering of the current. For a (5- function barrier, 
the filtering effect requires a non-magnetic scattering po- 
tential at the interface. This is not the case for a finite 
magnetic barrier due to the spin splitting of the Fermi 
surface. Specifically, minority-spin quasiparticle propa- 
gation is suppressed in the FM region above the critical 
injection angle 6>crit defined in Eq. (18), which is equiv- 
alent to a lowered transmissivity of the minority spins 
compared to the majority spins, giving rise to spin filter- 
ing. 

The z component of the spin current shows the a and (j) 
dependence expected for a spin- flipping mechanism. ^^'^^ 
This occurs due to the spin-dependent phase shift ac- 
quired by a triplet Cooper pair when it undergoes a spin- 
flip during the tunneling across the barrier. This mecha- 
nism requires a component of the magnetic moment par- 
allel to the d- vectors of both TSCs. In the perturbation 
theory for a (5-barrier junction the z component of the 
spin current satisfies^^ 



Tr2;,pert 



(X (Tsf cos (j) + 7i^sf ) sin 2a , 



(32) 



where Tgf and R^i are tunneling matrix elements of spin- 
flip tunneling and spin-flip reflection and 7 is an orbital- 
dependent factor connected to the phase shift experi- 
enced by specularly reflected quasiparticles. For the Py 
junction, there is no phase shift upon reflection and it is 



found that 7 = 1. In a junction, reflected Cooper pairs 
experience a tt phase shift leading to a factor of 7 = —1 
and in the more complicated case of a + ipy junction 
the phase shift of the Cooper pairs depends on their in- 
cident angle as tt — 2 diTctdin{ky/kz) so that — 1 < 7 < 0. 

We obtain the (/)-independent contribution /| q by per- 
forming a Fourier transformation of /| (</>)• Figure 10 
shows that /| q oscillates with increasing barrier width 
L for all considered gap symmetries. These oscillations 
are due to interference between quasiparticles reflected 
off the two interfaces. The sign difference between /| q 
for the py-waye and p^-wave cases agrees with previous 
results;^^'^^'^^ in particular, we find that the sign of the 
7 factor in Eq. (32) for the finite width barrier is the 
same as for the (5-function barrier. The {pz + ipy)-w8ive 
case shown in Fig. 10 is intermediate between these other 
cases, and we find that 7 can be positive or negative de- 
pending on A and L, in contrast to the (5- function results 
for which the sign was found to be always negative. 
Specifically, in Fig. 10 we see that for {pz + ipy)-wdive 
symmetry /| q is always positive at A = 0.3, at A = 0.5 it 
is negative for some values of L, and it is mostly negative 
at A = 0.8. 

In analogy to the procedure for the charge current, we 
demonstrate the appropriateness of Eq. (32) to describe 
the relevant physics, by showing that the oscillations of 
/| Q are due to oscillations of the spin- flip reflectivity T^sf • 
The parameter i?sf in Eq. (32) is expected to be roughly 
proportional to the reflectivity T^sf-^^ Since increasing 
the magnetization of the barrier makes it more likely for 
particles to be spin- flip reflected, i.e., T^sf should be en- 
hanced, we see that /| q increases with A for all gap sym- 
metries, see Fig. 10. We have also compared /| q to the 
weighted and ^-integrated spin-flip reflectivity 



donsf{o)Al 



(33) 



for a finite-barrier normal met al-FM- normal metal junc- 
tion. The weighting is carried out in the same way as 
for the charge current except for the Pz + ipy-waxe case, 
where it is necessary to account for the ^-dependent 
phase shift of quasiparticles reflected at the interface. 
Thus we weight the normal metal-FM-normal metal re- 
flectivity by Aq = sin^ — cos^ 0. Fig. 11 shows that the 
variation of /| q is in very good agreement with i^int, es- 
pecially for smaller values of L and A. We also find that 
the sign change in /| q for the {pz-\-ipy)-wdive case reflects 
an important change in the nature of the spin-flip reflec- 
tion:^^'^^ the positive /| q at small A indicates Py-like 
behaviour, where tunneling of quasiparticles with high 
incident angles ^ ~ 7r/2 is dominant. At higher values of 
A, however, the transmission of quasiparticles with lower 
incident angles 6 ^ is favored, leading to a p^-like be- 
havior and hence negative /| q . We hence conclude that 
when the L and A dependence of the transmissivity and 
reflectivity are taken into account, Eq. (32) qualitatively 
captures the observed behavior of the spin current. 
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FIG. 11: 0- independent part of the z component of the spin current (black sohd hne, left axis) and integrate spin- flip reflectivity 
Rint (red dashed line, right axis) as functions of the barrier width L for all considered gap symmetries. For p^^-wave the current 
is multiplied by —1 for better comparability. We set T — 0.4 Tc and A = 0.5 in all panels. 



IV. SUMMARY AND DISCUSSION 

In this paper we have studied the charge and spin cur- 
rents in a TSC-FM-TSC junction. For the charge cur- 
rent, we have found that transitions from a 0- to a tt- 
state are induced by rotating the magnetization of the 
FM barrier in the plane parallel to the d- vectors of both 
TSCs or by varying its magnitude. These effects have 
already been predicted for a J-function barrier. In ad- 
dition, we have found that changing the barrier width 
can also cause a O-tt transition. In all cases, O-tt tran- 
sitions are accompanied by sharp minima of the critical 
current, which can thus serve as diagnostics. For the case 
of {pz + ipy)-wdive symmetry, our results for the critical 
charge current agree with Rahnavard et al.^'^ However, 
we have gone beyond their work by considering two ad- 
ditional gap symmetries, by obtaining O-tt phase maps, 
and by extending the range of the exchange splitting. 
Furthermore, we have explained the origin of the charge 
current in terms of the variation of the spin-flip and spin- 
preserving transmissivities: for all considered gap sym- 
metries and barrier parameters, the latter dominates over 
the former in a 0-state, whereas the opposite is true for a 
TT-state. This result is consistent with lowest-order per- 
turbation theory. 

The spin currents in the finite-barrier case can be clas- 



sified in the same manner as for the J- function barrier, 
i.e., we find spin currents due to spin filtering and to spin 
flipping. Unlike for the (5-function barrier, however, the 
spin-filtering effect can occur in the absence of a finite 
charge scattering potential because of the spin splitting 
of the Fermi surface in the FM barrier. For the spin cur- 
rent in the {pz + ipy)-wdive case we find very different 
polarizations and phase dependences compared to Rah- 
navard et al.^'^ Specifically, the authors of Ref. 32 found 
a spin current polarization parallel to the d- vector, which 
is inconsistent with the fact that the spin of the Cooper 
pairs is always perpendicular to the d- vector. The ori- 
gin of this discrepancy might be that, in contrast to our 
work, in Ref. 32 the spin current was evaluated in the 
FM region, where the spin current is not a conserved 
quantity and is position dependent. 

Overall, our results show that approximating the bar- 
rier as a (5-function gives a qualitatively correct account 
of the Josephson currents in a TSC-FM-TSC junction. 
This demonstrates the robustness and, since a metallic 
ferromagnetic barrier is realizable in a Josephson junc- 
tion, the experimental relevance of the previously pre- 
dicted phenomena. ^^'^^'^^ In particular, the sign reversal 
of the current upon varying the orientation of the mag- 
netic moment is found to be a robust test of triplet su- 
perconductivity. The inclusion of a finite barrier width 
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is nevertheless crucial for a quantitatively realistic de- 
scription, as the barrier width can by itself control O-tt 
transitions of the junction. 
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